Equilibrium spin currents: Non-Abelian gauge invariance and color diamagnetism in 

condensed matter 
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The spin-orbit (SO) interaction in condensed matter can be described in terms of a non-Abelian 
potential known in high-energy physics as a color field. I show that a magnetic component of 
this color field inevitably generates diamagnetic color currents which are just the equilibrium spin 
currents discussed in a condensed matter context. These dissipationless spin currents thus represent 
a universal property of systems with SO interaction. In semiconductors with linear SO coupling the 
spin currents are related to the effective non-Abelian field via Yang-Mills magnetostatics equation. 
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SO interaction is considered as an important ingredi- 
ent of spintronics [HQ as it allows to control spin de- 
grees of freedom by electric means. Despite an increas- 
ing interest and growing number of publications on spin 
dynamics and spin currents in systems with SO interac- 
tion, a few basic questions remain unresolved up to now. 
In fact, even the very definition of spin currents is still 
debated 0, El, 0] • The reason for the controversies is that 
the spin is not conserved in a usual sense if SO interac- 
tion is present. The time derivative of the spin density 
s(r, t) can not be represented in form of divergence of 
a current, but always contains an extra term - the spin 
torque. Hence it appears that any redefinition of the cur- 
rent can be compensated by correcting the torque in a 
way that preserves dts. The problem of ambiguity of the 
spin current was sharpened by Rashba who noticed the 
presence of spin currents in a thermodynamically equi- 
librium 2D electron gas with Rashba SO interactions [5( . 
The physical reality of these dissipationless currents has 
been questioned as they do not accompanied by any spin 
accumulation, and their very appearance has been at- 
tributed to the ambiguity of the spin current concept. 

Recently @, @] Sonin proposed a way to detect the 
equilibrium spin flows in a " Rashba medium" , and ar- 
gued that there is nothing miraculous or unique in the 
presence of such equilibrium currents. They are known 
to exist in many different physical situations, e. g., spin 
supercurrents in a superfluid He 3 or Meissner currents 
in superconductors (see. Ref. @ and references therein). 
I would like to add another simple example that, as we 
will see, is much more relevant - the diamagnetic currents 
responsible for the Landau diamagnetism in metals. 

This paper is aimed at completely removing a flavor 
of ambiguity and mystery in the questions of equilib- 
rium spin flows and the definition of spin currents. The 
problem is resolved by making a link to non-Abelian 
gauge theories, and exploiting SU(2) gauge invariance 
of a many-body theory with SO interactions. The in- 
terpretation of SO interaction together with an external 



magnetic field as components A^ of a non-Abelian four- 
potential is known for many years [8| , and it is becoming 
more and more popular nowadays |9j|. However, the full 
power of non-Abelian gauge invariance in SO context is 
by far not explored [10]. In the present work I use this 
analogy to discover a simple physics behind equilibrium 
spin currents. We will see that these are nothing but dia- 
magnetic color currents that appear as a response to an 
effective Yang-Mills magnetic field produced by SO in- 
teraction. Due to the gauge invariance they are nonzero 
only if the field strength tFij is nonvanishing, which is the 
case in most real situations. This simple physical picture 
clearly demonstrates the universality of equilibrium spin 
currents in matter. They should be generically present 
in almost any system, like molecules or solids, if the SO 
coupling is non-negligible. Interestingly, in semiconduc- 
tors with linear SO coupling the spin current is related 
to the non-Abelian field by a Yang-Mills magnetostatics 
equation, which makes one more unexpected connection 
between the condensed matter and high energy physics. 
The present results also apply to ultracold atomic gases 
where a background non-Abelian field can be generated 



optically [ll|,|l2|. 



In general the many-body Hamiltonian with first spin- 
dependent relativistic corrections can be represented in 
a form M 



H = 



+ J7*t$ + i J dr'V rr ,&(r)&(r')y(r')V(r)} (1) 

where ^ = (tppt/j}) is a two-component fermionic field, 
V rr > is the interparticle interaction, and U is the ex- 
ternal scalar potential The components A^ (p = 
0,x,y,z) of the gauge field are 2x2 matrices of the form 
A-n = AuT a , where r a = a a /2 are the generators of 
SU(2) (spin-1/2 operators) with the following algebra 
[r a , T fc ] = ie ahc r c , tr{r a ,T 6 } = S ab 



The non-Abelian potential «4. M captures all spin effects 
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if one makes the following identifications 



Aq 



eh 



-B a T a 



At 



eh 



mc 



mc- 



(2) 



where B a are the components of the external magnetic 
field, and Ej is the electric field produced, for exam- 
ple, by nuclei in molecules or solids. The Hamiltonian of 
Eq. (fTJ) also covers popular 2D models of semiconductors 
with linear SO interaction of Rashba and/or Dresselhaus 
form |l5|. In this case the time component of SU(2) po- 
tential is still given by Aq — g[i^,B a , but with an appro- 
priate g-factor, while the spatial components are defined 
as follows, A z = 0, 

A x =2m{f3T x -o>t v ), Ay = 2m(ar x - (3t v ), (3) 

where a and (3 are the Rashba and Dresselhaus SO cou- 
pling constants respectively. 

The beauty of the representation |T|) is that the corre- 
sponding action 



S[V,A, 1 ] = dtidr&idt* - H) 



(4) 



is invariant with respect to local non-Abelian gauge 
transformations 



* i — ► U^>, 



Ap, i * UAJA- 1 - i^d^U- 1 (5) 



where U = e lSa<yr,t ^ ra is an arbitrary SU(2) matrix. The 
above gauge invariance immediately implies covariant 
conservation of a color current, J ^ = J^T a , with com- 
ponents J£ = 8S/5A": 



A Jo + DiJi = 



(6) 



where D^- — ■ — i[A lll •] is a covariant derivative, and 
SS 



(7) 



8Af 2m l v * ' 1 Am w 

Explicitly the covariant conservation law of Eq. (JSJ) reads 

d t J a + e abc A b J C + diJ? + e abc A\j^ = 0. 

Apparently the second and the fourth terms in this equa- 
tion violate conservation of the spin Jg = s a . The sec- 
ond term causes the spin precession in the U(l) magnetic 
field „4o- The fourth term is the "internal torque" due 
to SO interaction. However, the variational definition 
of the spin four-current, Eqs. (0 and ([8]), based on the 
gauge invariance, leaves no room for an ambiguity. As 
soon as we identify the zeroth component, SS/SAq, with 
the spin density, we are forced to accept that the spa- 
tial part, 5S/SAj, is the spin current. Jf is coupled to 
A1 in exactly the same fashion as Jg is coupled to Aq, 



which is absolutely analogous to the familiar case of the 
charge four-current coupled to U(l) gauge field. One can 
also show that J", Eq. ((HI, is a proper dissipative current 
conjugated to an effective SU(2) electric field J-^ [3]. It 
is worth noting that Eq. (jHJ) coincides with the "natural" 
definition of the spin current @, 0> 01 • 

Armed with the gauge invariant Hamiltonian and the 
variational definition of J" we are ready to approach 
the problem of equilibrium spin currents. SO interac- 
tion enters the Hamiltonian as an effective background 
non-Abelian field. If a magnetic part of this color field 
is nonzero one naturally expects an orbital response in 
a form of color diamagnetic currents. These currents, 
if exist, are given by the derivative of the energy (ther- 
modynamic potential) S[-4°] = (H) with respect to Af. 
Since the energy is gauge invariant it can depend on Af 
only via invariants composed of the field strength 



'ij — &i iA-j Oj >Ai I \>Ai , >A-j ] 



(9) 



A particular form of invariants is determined by the sym- 
metry of a particular system. 

For the sake of clarity I consider explicitly the case of 
semiconductors with linear SO coupling of the Rashba- 
Dresselhaus type. I also assume that the external scalar 
potential and the usual magnetic field are zero, i. e. in 
Eq. (fT]) U — and Aq = 0. Since in the absence of Af the 
system is rotationally invariant, the first SO correction 
to the energy must be proportional to tr{TijTij), i. e. 



\ I 'WW. 



(10) 



where A is a constant (on dimensional grounds A ~ 
Pp~ 2 /m, where pf is the Fermi momentum and D > 1 
is the dimension of space). Calculation of the current, 
Jf = -SE so /SAf, yields 



Jj — XDi^F ij 



XidiFij-iiAi,^]), 



11) 



which is exactly of the form of Yang- Mills magnetostat- 
ics equation. Physically the result is very similar to the 
case of U(l) magnetic field: an external field produces 
diamagnetic currents aimed at compensating that field. 
There is, however, an essential difference in the spa- 
tial distribution of diamagnetic currents. In U(l) case 
the currents in the bulk vanish when the magnetic field 
approaches a constant. Only the integral defining the 
induced magnetic moment remains finite. In contrast, 
in the non-Abelian case bulk diamagnetic currents ex- 
ist even for a constant in space field. The reason is the 
commutator in the right hand side of Eq. (jlljl . A simi- 
lar commutator in Eq. ^ gives a nonzero magnetic field 
even for a space-independent vector potential. Thus, in 
the case of a homogeneous field (space-independent SO 
coupling constants), the bulk spin current is given by 



-X[Ai, [Ai,Aj}} 



(12) 
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Equation (fT2|) shows that in the homogeneous system the 
spin current is proportional to the third power of non- 
Abelian potential. This naturally explains why Rashba's 
equilibrium spin current is proportional to a 3 Q . 

The above phenomenology can be confirmed by direct 
microscopic calculations of the spin current for an ex- 
actly solvable model of noninteracting particles. First I 
consider homogeneous Ai- In this case the spin current, 
Eq. ((8]), is given by the expression 



J? = T 



4m 



Af 



(13) 



where cj is the fermionic Matsubara frequency, and n is 
the density of particles. The one particle Green's func- 
tion G(u!,p) is defined as follows 



G(w,p) = 
Z+(P) 



^ + M-^(Pi-^r fe ) 2 
Z-(P) 



fi-E + (p) iuj + (i - E- (p) ' 



where 



^( Pj A b )( Pk A b 



(14) 



(15) 



P 2 , AfAf , 1 



^VfeA^M) ( 16 ) 



2m 8m 2m 
After summation over uj Eq. (TIB"]) takes the form 

J? = ^-E [n F {E-)-n F {E +) ] 

) 771 ' 



2m 



J2[nF(E-) + n F (E + )}, 



(17) 



where n F (E) is the Fermi distribution function. In the 
second term in Eq. (fT7|) the density of particles is rep- 
resented as n = J2 p [ n F(E-) + n F (E + )]. The rest of 
calculations is straightforward. Assuming as usual that 
SO coupling is weak, we expand the distribution func- 



tions n F (E±) in terms of J (pjA b )(pkA b k ), and keep the 

first nonvanishing term in Eq. (|17[). This term is, as ex- 
pected, of the third order in Aj. The final result at zero 
temperature is the following 



J? 



N F 
24m 2 



( A a A b A b - A a A b A b \ 



(18) 



where Np is the density of states at the Fermi level. It is 
easy to see that the expression in the brackets in Eq. (fT8|) 
is exactly the double commutator entering the right hand 
side of Eq. (fl"2"]) . Hence the phenomenological coefficient 
in the SO energy is A = N F /24m 2 . Thus the direct 
calculations indeed confirm a diamagnetic nature of the 



equilibrium spin currents. The spin current in the 2D 
Rashba-Dresselhaus model is obtained by inserting Af 
of Eq. © and N F = m/2-n: into Eq. (□ 



(19) 
(20) 



Setting f3 = we exactly recover the result by Rashba 
[5j. It is very interesting to realize that the formula for 
the spin current obtained in Ref. is a hidden form of 
a covariant curl of the non-Abelian magnetic field! The 
spin current of Eqs. (p~9|) . (|20|) is zero at a = ±/3. The 
reason is that the color magnetic field T^ y = 4m 2 (a 2 — 
(3 2 ) vanishes at these special values of SO constants. In 
the absence of magnetic field there are no diamagnetic 
currents. 

Now I will show that the gauge invariance allows to 
significantly simplify practical calculations of the spin 
current using powerful techniques of the linear response 
theory, in spite of an obvious nonlinearity in Af . Let us 
first assume that the non-Abelian potential Af (r) is weak 
(here we allow for a general inhomogeneity of the poten- 
tial). The standard linear response theory [ItJ yields 



•TO 



J dr X ? 6 (r,r')^(r'), 



(21) 



where the response function is defined by the Kubo for- 
mula 



X^(r,r') = «Jf(r);j}(r')))-5(r-r')^ ( 5 afc £. (22) 



The operator Jf in Eq. (j22|) is given by the first term in 
the right hand side of Eq. ([5]) . 

Due to the gauge invariance the vector potential Af 
can enter Eq. (f2Tj) only via the linearized field strength, 
diA°j — djAf, and, possibly, its spatial derivatives, <9,. By 
the same token, to recover the formula valid to the first 
order in the full nonlinear field strength, of Eq. (|11[) . 
we simply have to make the following replacements 



(cUf - d 3 Af) 



0, 



D, 



(23) 



All other changes are forbidden by the gauge invariance. 

To demonstrate how this prescription works in practice 
I consider again the example of a homogeneous (in the 
absence of Af(r)) and noninteracting electron gas with 
linear SO interaction. In this case Eq. (|22[) for the spin 
current response function in the momentum representa- 
tion simplifies as follows xfH^l) — ^ afc X^(q)i 



x|-(q) = 



n f (g p _a) - w F (£ p+ a) 
2m 2 (£ p+a -£ p _ f ) 



4m 



(24) 



where £ p = Equation ([24)1 coincides, up to a numer- 
ical factor, with a static charge current response function 
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Xij (q) 17] . It the standard theory the function Xij (q) de- 
termines the Landau diamagnetic response to an external 
magnetic field, which provides us with another confirma- 
tion of the diamagnetic nature of the equilibrium spin 
currents. To calculate the spin current response function 



we follow the usual route [17|. In the limit of q <C Pf, 



which means that SO constants are slowly changing on 
the scale of pp 1 , x«(q) takes the form 



Xij (q) rt A 2 

J 24m z 



5 a \q iq] -q 2 5 l3 ). 



Transforming this equation from q- to r-space, and in- 
serting the result into Eq. (j2Tj) we get for the spin current 



7° 



N F 
24m 2 



(25) 



The final step is a substitution of Eq. (j2li|) , which yields 
the final gauge covariant expression valid to the first or- 
der in the total non-Abelian magnetic field Tij, Eq. 



3 24m 2 



D T a 



Thus the standard linear response theory supplemented 
with the substitution of Eq. (j23|) exactly recovers the 
Yang-Mills form of the spin current, Eq. (fTTj) . with the 
correct coefficient A = Nf /24m 2 . An obvious advantage 
of this way is that it allows to straightforwardly include 
the effects of interaction and external inhomogencitics 
(for example, impurity scattering). The formalism is also 
easily transferable to nonequilibrium situations. 

Apparently the general conclusion about the nature of 
equilibrium spin currents is not restricted to the simple 
Rashba-Dresselhaus model. Such diamagnetic currents 
necessarily present in any system with SO interaction, 
provided the effective color magnetic field of Eq. ([9]) is 
nonzero. As vanishing !Fij is actually an exception (like 
in the case a — ±(3) it should be possible to find the 
spin currents in many molecules and solids using available 
codes for ab initio electronic structure calculations. 

In conclusion I identified the equilibrium spin currents 
with diamagnetic color currents in the presence of a non- 
Abelian field generated by SO coupling. If the particles 
have a color charge coupled to physical Yang-Mills fields, 
like in quark-gluon plasma, the color currents would pro- 
duce a back reaction field to compensate the external 
one. The absence of such back reaction does not make 
the equilibrium spin currents less " real" . They are as 
real as pseudo-diamagnetic currents in a rotating refer- 
ence frame. One of the main outcomes of this work is a 
conclusion about universality of equilibrium spin currents 
that should exist in most real systems. The universality 
makes it especially intriguing to observe and, possibly, 
to control such currents experimentally. In this respect 
a proposal by Sonin Q > connecting spin currents to me- 
chanical deformations, looks especially interesting. The 



present results show that it is not necessarily to do exper- 
iments with a "Rashba medium". Any technologically 
convenient material with strong SO interaction should 
demonstrate the same effect. Possibly the color diamag- 
netic currents can be also observed in trapped atomic 
gases subjected to aproper configuration of artificial non- 



Abelian fields [111, [12 1 
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